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Abstract. Recent in situ and satellite measurements sug- 
gest a contribution of ~5 pptv to stratospheric inorganic 
bromine from short-lived bromocarbons. We conduct a mod- 
eling study of the two most important short-lived bromocar- 
bons, bromoform (CHBo) and dibromomethane (CH^B^), 
with the Goddard Earth Observing System Chemistry Cli- 
mate Model (GEOS CCM) to account for this missing 
stratospheric bromine. We derive a “top-down” emission 
estimate of CHBrs and CHaBra using airborne measure- 
ments in the Pacific and North American troposphere and 
lower stratosphere obtained during previous NASA aircraft 
campaigns. Our emission estimate suggests that to repro- 
duce the observed concentrations in the free troposphere, a 
global oceanic emission of 425GgBryr _l for CHBrj and 
57 GgBryr~' for CHiBri is needed, with 60% of emissions 
from open ocean and 40% from coastal regions. Although 
our simple emission scheme assumes no seasonal variations, 
the model reproduces the observed seasonal variations of the 
short-lived bromocarbons with high concentrations in win- 
ter and low concentrations in summer. This indicates that the 
seasonality of short-lived bromocarbons is largely due to sea- 
sonality in their chemical loss and transport. The inclusion 
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of CHBri and CH 2 Br 2 contributes ~5 pptv bromine through- 
out the stratosphere. Both the source gases and inorganic 
bromine produced from source gas degradation (Br 3 ’ SLS ) in 
the troposphere are transported into the stratosphere, and 
are equally important. Inorganic bromine accounts for half 
(2 .5 pptv) of the bromine from the inclusion of CHBr 3 and 
CH 2 Br 2 near the tropical tropopause and its contribution 
rapidly increases to ~ 100% as altitude increases. More than 
85% of the wet scavenging of Br 3SLS occurs in large-scale 
precipitation below 500 hPa. Our sensitivity study with wet 

scavenging in convective updrafts switched off suggests that 
Br vsLS 

in the stratosphere is not sensitive to convection. 
Convective scavenging only accounts for ~0.2 pptv (4%) dif- 
ference in inorganic bromine delivered to the stratosphere. 


1 Introduction 

Oceanic emission of very short-lived substances (VSLS) is 
thought to contribute significantly to reactive bromine in 
the stratosphere in addition to long-lived halons and methyl 
bromide (Kurylo and Rodriguez, 1999). VSLS are not ac- 
counted for in most chemistry climate models. In the strato- 
sphere, inorganic bromine produced from VSLS (Br 3SLS ) 
contributes to catalytic destruction of ozone (e.g., McElroy et 
al., 1986; Solomon et al„ 1995; Garcia and Solomon, 1994; 
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Sturges et al., 2000). Br y SLS can also have a significant im- 
pact on tropospheric ozone (von Glasow et al., 2004; Yang et 
al., 2005). 

A key question is the extent to which VSLS contribute 
to inorganic bromine (Br y ) in the stratosphere. The dif- 
ference between ground-based observations of column BrO 
(Sinnhuber et al., 2002, and references therein) and the frac- 
tion of the BrO column which can be accounted for based 
on long-lived halons and methyl bromide suggests a contri- 
bution of ~5 pptv to stratospheric inorganic bromine from 
VSLS (Sinnhuber et al., 2002). This is consistent with recent 
balloon-borne measurements which estimated that Br^ SLS is 
~5.2pptv (Dorf et al., 2008). Estimates of Br^ SLS implied 
from satellite BrO measurements are more variable, ranging 
from ~3 pptv (Sinnhuber et al., 2005; Livesey et al., 2006) 
to ~8pptv (Sioris et al., 2006). The estimated contribution 
of VSLS to Br y in the stratosphere from modeling studies 
are in general significantly lower than that suggested by in 
situ and satellite observations. Earlier studies by Dvortsov 
et al. (1999) and Nielsen and Douglass (2001) calculated a 
maximum of 1.8 pptv and 1.1 pptv Br v , respectively, in the 
stratosphere from bromoform (CHBrj). More recent model- 
ing studies by Kerkweg et al. (2008), Gettelman et al. (2009), 
Aschmann et al. (2009) and Hossaini et al. (2010) suggested 
that CHBr 3 and CH 2 Br 2 contribute significant amounts of 
reactive bromine to the lower stratosphere, with the calcu- 
lated contribution varying between ~2-3 pptv among indi- 
vidual studies. One exception is the modeling study by 
Warwick et al. (2006) which presented a detailed emission- 
based modeling analysis of all five major VSL oceanic bro- 
mocarbons, including CHBr 3 , dibromomethane (CHiBn), 
bromodichloromethane (CHBrCL), dibromochloromethane 
(CHBnCl), and bromochloromethane (Cb^BrCl). They 
concluded that very short-lived (VSL) oceanic bromocarbons 
contribute a maximum of 6-7 pptv Br y , peaking over the 
equator at ~ 100 hPa. 

A large part of the uncertainty in the contribution of VSLS 
to Br y in the stratosphere arises from the highly-variable 
mixing ratios of VSL bromocarbons in the marine bound- 
ary layer, ranging from 0.2 pptv to >100 pptv (Quack and 
Wallace, 2003). Therefore, the amount of VSL bromocar- 
bons that enters the free troposphere and stratosphere and its 
subsequent inorganic bromine may vary depending on the 
location and timing of emissions. The current estimate of 
the oceanic emissions of VSL bromocarbons is highly un- 
certain. Bottom-up emission estimates based on atmospheric 
and oceanic surface water observations indicate that sea-to- 
air flux of CHBr 3 is between 210-820Gg Bryr” 1 (Carpen- 
ter and Liss, 2000; Quack and Wallace, 2003; Yokouchi et 
al., 2005; Butler et al., 2007) and that for CbbBrj ranges 
between 61-280 Gg Bryr” 1 (Yokouchi et al., 2005; Butler 
et al., 2007). Modeling studies that derive top-down emis- 
sion estimates by reproducing background atmospheric con- 
centrations suggest that global emission of CHBn ranges 
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between 1 90-570 Gg Bryr” 1 (e.g. Dvortsov et al., 1999; 
Kurylo and Rodriguez, 1999; Nielsen and Douglass. 2001; 
Warwick et al., 2006) and that of CFLBrs is between 43- 
104GgBryr”' (e.g. Dvortsov et al., 1999; Warwick et al., 
2006). 

Another important contributor to the uncertainty is re- 
lated to the transport and dehydration process in the tropi- 
cal tropopause layer (TTL). Convective lofting is suggested 
to be the most important pathway of air entering the strato- 
sphere (e.g. Sinnhuber and Folkins, 2006; Fueglistaler et al., 
2009). Air masses entering the stratosphere must be dehy- 
drated and soluble substances, e.g. Br y , carried within these 
air masses are subject to wet scavenging during troposphere- 
to-stratosphere transport. The process of dehydration is not 
well understood. Past studies suggest that dehydration can 
occur slowly by falling ice during slow ascent through the 
TTL (Holton and Gettleman, 2001 ; Fueglistaler et al., 2005) 
or more rapidly in deep convective overshooting (Sherwood 
and Dessler, 2000). As a result, the rate of wet scaveng- 
ing of Br y , which depends on the mechanism of dehydration, 
is highly uncertain. Model sensitivity studies suggest that 
the contribution of CHBn to stratospheric bromine range 
from 0.5-1. 6 pptv when Br y is removed completely in con- 
vective updraft to ~3 pptv assuming no washout (Sinnhuber 
and Folkins, 2005; Aschmann et al., 2009). 

In this paper, we use the GEOS Chemistry Climate Model 
(GEOS CCM) together with aircraft measurements from pre- 
vious NASA field missions to derive a top-down emission es- 
timate of CHBr 3 and CH 2 Br 2 , the two most important VSL 
bromocarbons. Together, CHBr 3 and CH 3 Br 2 account for 
>80% VSL organic bromine in the marine boundary layer 
and free troposphere (WMO, 2007). We present a quanti- 
tative estimate of their contribution to reactive bromine in 
the stratosphere with our optimized emission. Compared to 
many previous studies, this study has two major improve- 
ments. First, we derive the top-down emission estimate 
and evaluate model results with an extensive set of aircraft 
measurements obtained during eight previous NASA field 
missions throughout the Pacific/North America region. The 
broad spatial and temporal coverage of these measurements 
lends more confidence in the robustness of our derived re- 
sults. Secondly, unlike many previous modeling studies 
that assume a constant wash-out lifetime of Br y (10-days in 
Dvortsov et al., 1999, Nielsen and Douglass, 2001 ; Warwick 
et al., 2006; Hossaini et al., 2010; 0 to indefinite days in 
Sinnhuber and Folkins, 2005) or a relative simple washout 
process (Aschmann et al., 2009), we implement a detailed 
wet deposition scheme that includes scavenging in convec- 
tive updrafts and rainout/washout in large-scale precipitation. 
This is to assure a more realistic representation of the impact 
of large-scale/convective transport on the highly soluble Br y . 
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2 Model 

2.1 Description 

Model simulations are conducted using the Goddard Space 
Flight Center (GSFC) GEOS CCM Version 2 (V2), an aug- 
mented version of GEOS CCM Version 1 ( V 1 ) as described 
in Pawson et al. (2008). The model has a spatial resolu- 
tion of 2° latitude by 2.5° longitude and 72 layers extending 
from the surface to 0.01 hPa. The model couples the GEOS- 
5 GCM (Reinecker et al., 2008) with an updated version of 
the stratospheric chemistry module described by Douglass 
and Kawa (1999). The photochemical scheme includes all 
important gas phase reactions for the stratosphere (Douglass 
and Kawa, 1999) and uses the chemical kinetics from JPL 
2006. Photolysis rates are calculated using the temperature 
dependent cross sections from JPL 2006. The model uses 
a flux-form semi-Lagrangian dynamical core (Lin, 2004). 
Moist processes in GEOS-5 are represented using a convec- 
tive parameterization and prognostic cloud scheme. Convec- 
tion is parameterized using the relaxed Arakawa Schubert 
(RAS) scheme developed by Moorthi and Suarez (1992) in 
which the atmosphere is relaxed towards equilibrium. Ott et 
al. (2010) provides a detailed analysis of the GEOS-5 moist 
physics scheme and its impact on tracer transport, describing 
the prognostic cloud scheme and comparing results with ob- 
servations. The scheme calculates large-scale ice and liquid 
condensate by assuming a probability distribution function 
for total water. Condensate is subsequently removed from the 
model domain by evaporation, auto-conversion of liquid con- 
densate, sedimentation of frozen condensate, and accretion 
of condensate by falling precipitation. The GEOS-5 GCM 
produces general precipitation patterns that correlate well 
with the Global Precipitation Climatology Project (GPCP) 
data (r =0.65). The stratospheric and tropospheric transport 
produced by the GEOS-CCM has been extensively evaluated 
in many previous studies (e.g. Eyring et al., 2006; Waugh et 
al.,2007; Pawson et al.,2008; Douglass et al., 2008; Liang et 
al., 2008). The GEOS CCM simulations agree well with ob- 
servations in many of the meteorological, transport-related, 
and chemical diagnostics. Of particular relevance to this 
study, Douglass et al. (2008) and Waugh et al. (2007) show 
that the GEOS CCM reproduces well the mean atmospheric 
circulation with realistic age-of-air, and therefore a realis- 
tic representation of the troposphere-to-stratosphere and the 
subsequent transport. 

2.2 Bromine chemistry in the GEOS CCM 

For this study, we modify the standard stratospheric pho- 
tochemistry scheme to include simple VSL bromocarbon 
chemistry. In this simple chemistry scheme, CHBrj and 
CFLBri are destroyed via photolysis and reaction with OH. 
Photolysis rates are calculated using the JPL 2006 cross 
sections. Rates of reaction with OH are calculated using 
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temperature-dependent expressions from JPL 2006. OH 
above the tropopause is calculated online in the stratospheric 
chemistry module. Below the tropopause, OH is relaxed to 
zonal-averaged 2-dimensional monthly mean fields archived 
during a full chemistry simulation from the GEOS-Chem 
chemistry transport model (Bey et al., 2001). GEOS-Chem 
includes detailed tropospheric 03 -NO x -NMHC chemistry 
and the simulated OH compares reasonably well with obser- 
vations (Bey et al., 2001). The mean tropospheric OH con- 
centration is 1 1 x 10 5 molecules cm~ 3 and yields an atmo- 
spheric lifetime of 5.1 years for methyl chloroform, within 
the uncertainty of observational estimates (Bey et al., 2001). 

Since >80% of the inorganic bromine produced by the 
degradation of CHBrj and CHiBr 2 (Br y SLS ) in the tropo- 
sphere is in the form of hydrobromic acid (HBr) and hy- 
pobromous acid (HOBr) and both are highly soluble (Yang 
et al., 2005), for simplicity, we track the degradation prod- 
ucts as a single tracer. Br y SLS is transported as an individual 
tracer and subject to wet and dry deposition, but does not in- 
teract with standard stratospheric chemistry. We adopt the 
wet and dry deposition scheme from the Goddard Chemistry 
Aerosol Radiation and Transport (GOCART) model (Chin 
et al., 2000a). Dry deposition includes gravitational set- 
tling with a uniform deposition velocity of 0.3 cm/s, a sim- 
plified value based on Yang et al. (2005). Wet deposition 
includes scavenging in rainout (in-cloud precipitation) and 
washout (below-cloud precipitation) in large-scale precipita- 
tion (Giorgi and Chameides, 1986) and in deep convective 
updrafts (Balkanski et al., 1993). The VSL bromocarbons 
are in general insoluble and not subject to wet scavenging. 
We assume high solubility for Br y . Washout and rainout of 
Br y by large-scale rain are computed as first-order processes 
using parameters from Giorgi and Chameides (1986) and the 
simulated precipitation rate. Convective scavenging occurs 
in the moist convection process with a scavenging efficiency 
of 1 km“' , equivalent to 100% removal. When evaporation 
occurs during large-scale and convective transport, a frac- 
tion of the dissolved Br y proportional to the evaporated wa- 
ter is released back to the atmosphere (Chin et al., 2000a). 
This wet scavenging scheme has been applied in many at- 
mospheric modeling studies of soluble aerosols, e.g. sulfate, 
sea salt, dust, and the simulated concentrations compare well 
with surface observations at many observation sites around 
the globe (e.g. Chin et al., 2000b, 2007; Ginoux et al.,2001). 

Figure 1 presents the calculated lifetime of CHBr 3 and 
CHjBra from the GEOS CCM bromocarbon simulation. Lo- 
cal lifetime of VSL bromocarbons varies significantly from 
~15 days for CHBr 3 and ~90 days for CH- Bri in the trop- 
ics to ~ 180 days (CHBr;) and ~540 days (CHiBra) in the 
poles, similar to the results in Dvortsov et al. (1999) and Hos- 
saini et al. (2010). Dvortsov et al. (1999) calculated that the 
tropical tropospheric lifetime of CHBri is of the order of 
2-3 weeks and that of CHiB^ is 2-3 months. Hossaini et 
al. (2010) showed that local lifetime of CHBr 3 ranges be- 
tween ~15 days (compared to ~13 days in this work) at 
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CHBr 3 Lifetime against photolysis (days) CHBr g Lifetime against OH (days) CHBr 3 Total lifetime (days) 



90°S 60°S 30°S 0 30°N 60°N 90°N 

Latitude 


90°S 60°S 30°S 0 30°N 60°N 90°N 

Latitude 


90°S 60°S 30°S 0 30°N 60°N 90°N 

Latitude 




Q°S 60°S 30°S 0 30°N 60°N 90°N 90°S 60°S 30°S 0 30°N 60°N 90°N 

Latitude Latitude 


CH 2 Br 2 Lifetime against photolysis (days) 

"T” 


Lifetime against OH (days) 


90°S 60°S 30°S 0 30°N 60°N 90°N 

Latitude 


Fig. 1. Model calculated local lifetime of CHBr3 and CfHB^ against photolysis, reaction with OH, and total lifetime. 


the tropical surface to ~25-30 days (30 days this work) in 
the TTL, and that of ranges from ~50 days (~70 

days this work) at the surface to ~520 days (~540 days 
this work) in the TTL. Though local lifetime of VSL bromo- 
carbons is highly variable, the global-integrated atmospheric 
lifetime is heavily weighted by loss in the tropical lower at- 
mosphere where emissions are high and degradation occurs 
rapidly. The global-integrated atmospheric lifetime (r) of 
CHBr-t is 20 days (annual average) in our simulation with 
photolysis being the dominant sink (r/ It ,=27 days, toh=81 
days), consistent with the lifetime calculated in Warwick et 
al. (2006) (varying between 15-37 days with the choice of 
emissions) and Kerkweg et al. (2008) (20 days). Degrada- 
tion of CH 2 Br 2 occurs predominantly via reaction with OH 
(t a „=8300 days, toh=143 days). The simulated atmospheric 
lifetime of O-DBm is 140 days, longer than the 100 days cal- 
culated in Kerkweg et al. (2008). Our simulated seasonal cy- 
cle of CHBr 3 at Hawaii (maximum in winter ~0.7 pptv and 
minimum in summer ' 0.3 pptv) matches well with obser- 
vations from Atlas and Ridley (1996). This seasonal cycle 
is generated by transport and chemistry in the lower tropo- 
sphere (Nielsen and Douglass, 2001), which implies that the 
model captures well chemical loss and tropospheric transport 
processes. 

The simulated atmospheric lifetime of Br y SLS against 
wet deposition is ~15 days with column-integrated rain- 
out/washout lifetime varying from ~10 days in the tropics 
to ~2 months at high latitudes. The calculated lifetime of 
Br y LS against dry deposition is greater than 20 years, sug- 
gesting surface dry deposition is of negligible importance 
compared with wet scavenging. 


3 Emission scenarios 

In this work, we derive a “top-down” emission scenario ( Sce- 
nario A) of CHBrs and CHiBn using observations from pre- 
vious aircraft campaigns as constraints. We conduct a 12- 
year GEOS CCM bromine simulation with emission Sce- 
nario A and examine the chemistry and transport of VSL 
bromocarbons and their degradation product, Br,y sl s , in the 
atmosphere. The simulation is driven by sea surface tem- 
peratures (SST) from 1990 to 2001 , but otherwise does not 
correspond to any specific year. We initialized the simula- 
tion with zero-concentrations for CHBo and O-LBn. For 
CHBr 3 (lifetime ~1 month) and CH 2 Br 2 (lifetime of ~4 
months) a 3 year simulation length has been used by most 
previous modeling studies and is generally thought to be ad- 
equate. However, transport in the stratosphere is slow and 
Br y has no significant loss process in the atmosphere other 
than removal by wet scavenging in the troposphere, thus a 3 
year simulation is not sufficient to quantitatively understand 
the contribution of VSL bromocarbons to reactive bromine 
in the stratosphere. Here we extend the simulation length 
to 12 years to investigate whether the duration of simulation 
affects the troposphere-to-stratosphere transport of Br y . 

We conduct two sensitivity simulations with two addi- 
tional emission scenarios (Scenario B and C) from Warwick 
et al. (2006), which were demonstrated to reproduce well 
the observed CHBr 3 concentrations during the NASA PEM- 
Tropics mission. For computational efficiency, the two sensi- 
tivity simulations are run for three years driven by SST from 
1990-1992. Results from year 3 from all three simulations 
are compared with each other to understand the impact of 
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Table 1 . Emission distribution of CHBrj and CH 2 Br 2 in Scenario A, B and C. 


Global emissions 
GgBryr" 1 

CHBr 3 (CH 2 Br 2 ) a 



Emission distribution 



Latitude 

Open ocean 
Percentage 

GgBryr la 

Latitude 

Coast 

Percentage 

GgBryr" 111 

Scenario A 

425 (57) 

5-90° S 

1.1% 

5(1) 

50-90° S 

1.1% 

5(1) 



10-50° S 

11.1% 

47(6) 

10-50° S 

11.1% 

47(6) 



10° S-10° N 

33.3% 

141 (20) 

10° S-10° N 

4.4% 

19(3) 



10-50° N 

13.3% 

57 (8) 

1 0-50° N 

22.2% 

94(13) 



50-90° N 

1.1% 

5(1) 

50-90° N 

1.1% 

5(1) 

Scenario B 

380(104) 

20° S-20° N 

75% 

285 (78) 

None 





20-50° S/N 

25% 

95 (26) 




Scenario C 

565 (57) 

20° S-20° N 

37.8% 

214(21) 

10° S-10° N 

49.6% 

280 (29) 



20-50° S/N 

6.3% 

71 (7) 





a Emissions of CH 2 Br 2 are in the parentheses, 

emission distribution on simulated bromocarbon concentra- 
tions. 

The details of each scenario are listed in Table 1 and emis- 
sion distributions are presented in Fig. 2. For simplicity we 
assume no seasonal variability in emissions in all three sce- 
narios. We refer to the three simulations as simulation A 
(with emission Scenario A), simulation B (with Scenario B), 
and simulation C (with Scenario C). 

3.1 Scenario A 

We use measurements from eight NASA aircraft missions be- 
tween 1996 and 2008 as observational. constraints. A detailed 
list of the aircraft missions is presented in Table 2 and the 
geographic distribution of all flight tracks is shown in Fig. 3. 
The composite of bromocarbon measurements in the tropo- 
sphere covers the entire Pacific and North America region. 
Measurements in the upper troposphere/lower stratosphere 
(UT/LS) are concentrated in the tropics. Whole air samples 
were collected by stainless steel canisters onboard the NASA 
DC-8, ER-2, and WB-57F aircrafts and later analyzed using 
gas chromatography with mass selective detection (GC/MS) 
(Schauffler et al., 1999; Blake et al., 2003). 

We derive the emission distribution of CHBr 3 following a 
similar approach as that introduced in Warwick et al. (2006) 
and use Scenario 5 of Warwick et al. (2006) (Scenario C be- 
low) as our baseline emission. We use 10 individual compo- 
nents to represent emissions from open ocean and coastlines 
near the equatorial tropics, and at the mid and high latitudes, 
respectively. Multiple sensitivity runs are conducted by vary- 
ing the magnitude of emission for each region and the lat- 
itudinal border of each region. The optimum emission esti- 
mate for each region is determined by simultaneously match- 
ing i) the observed concentration in the middle troposphere 
and ii) the observed vertical gradient in the corresponding 





(a) Scenario A: 425 

90 N 


60 N 
30°N 
0 ° 
30°S 
60°S 


90 S 
180° 


120°W 60°W 


180° 120°W 60°W 


180° 120°W 60°W 


60°E 120°E 


(10 kg/m'/s) 


Fig. 2. Global emission distribution of CHBr 3 for (a) Scenario A, 
(b) Scenario B, and (c) Scenario C. 
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Table 2. A summary of aircraft observations used in this study. 


Missions 

Time 

Area 

References 

PEM-Tropics 1 

PEM-Tropics A 
PEM-Tropics B 

August -October 1996 
M arch-April 1999 

75° S-50° N, 150° E-100° W 
40° S-40° N, 140° E-70° W 

Hoeilet al. (1999) 
Raper et al. (2001) 

TRACE-P 2 


February- April 2001 

10-40° N, 1 40° E-80° W 

Jacob et al. (2003) 

INTEX 3 

1NTEX-A 

INTEX-B 

July-August 2004 
March-May 2006 

10-45 °N. 70-130° W 
15-70° N, 170° E-90° W 

Singh et al. (2006) 
Singh et al. (2009) 

TC4 4 


July-August 2007 

10° S-40°N, 60-130° W 

Toon et al. (2010) 

ARCTAS 5 

ARCTAS-A 

ARCTAS-B 

April 2008 
June-July 2008 

30-90° N. 20° E-50° W 
30-90° N, 20-140° W 

Jacob et al. (2009) 
Jacob et al. (2009) 

STRAT 6 


January-December 1996 

5° S-65°N, 80-160° W 

Schauffleret al. (1999) 

Pre-AVE 7 


January-February 2004 

10° S-40° N, 80-100° W 


AVE 8 

HAVE-2 

CRAVE 

June 2005 

January-February 2006 

15-50° N, 70-100° W 
0-40° N, 80-100° W 

Kroon et al. (2008) 
Kroon et al. (2008) 


1 The Pacific Exploratory Missions - Tropics (PEM-Tropics) mission. 

2 The Transport and Chemical Evolution over the Pacific (TRACE-P) mission. 

3 The Intercontinental Chemical Transport Experiment (INTEX) mission. 

4 The Tropical Composition, Cloud and Climate Coupling (TC4) mission. 

5 The Arctic Research of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS) mission. 

5 The Stratospheric Tracers of Atmospheric Transport (STRAT) mission. 

7 The Pre-Aura Validation Experiment (Pre-AVE). 

8 The Aura Validation Experiment (AVE). Data from two AVE campaigns are used in this study: i) 2005 June Houston AVE (HAVE-2) 
campaign, ii) 2006 January Costa Rica AVE (CRAVE) campaign. 



60°E 120°E 180° 120°W 60°W 0° 


Fig. 3. Flight tracks for all aircraft missions used in this study (Ta- 
ble 2). 


region. Measured concentrations of CHBr 3 and CH 2 Br 2 in 
the coastal marine boundary layer are highly correlated (Yok- 
ouchi et al.. 2005), suggesting they are produced by sim- 
ilar marine macroalgae sources and therefore have similar 
source distributions (Carpenter et al., 1999). Atmospheric 
measurements from the aircraft campaigns used in this study 
also show consistent high correlations (correlation coeffi- 
cients range between 0.59 and 0.96 for individual missions) 
between the measured concentrations of CFIBr? and CH 2 Br 2 
below 1 km altitude. Therefore, we employ the same distri- 


bution for CH 2 Br 2 as that derived for CHBr 3 . The global 
emission magnitude for CH 2 Br 2 is deduced by matching the 
simulated concentrations with background atmospheric ob- 
servations. To avoid redundancy, we present here only the 
optimum emission estimate (Scenario A). 

In our optimum emission estimate, we separate global 
emissions into 10 individual components representing emis- 
sions from open ocean and coastlines at 10°S-T0°N, 10- 
50° S/° N, and 50-90° S/° N, respectively (Fig. 2 and Ta- 
ble 1). Our best estimate yields a global emission of 
425 GgBryr~* forCHBr 3 , with 255GgBryr~' (60%) and 
1 70 Gg Bryr - 1 (40%) over the open ocean and along the 
coastlines, respectively. Our estimate of CHBr 3 emission 
from coastal regions is similar to the bottom-up coastal es- 
timate of 200 Gg Br yr 1 from Carpenter et al. (2009). The 
emission is latitude-dependent, with 160GgBryr~~ 1 (~38%) 
concentrated in a narrow tropical band between 10° S and 
10° N,95 Gg Bryr -1 (22%) and ffiOGgBryr- 1 (35%) in the 
Southern and Northern sub-topics and mid-latitudes between 
10° and 50° north and south, respectively. A small amount, 
lOGg Bryr -1 , is emitted at the Northern and Southern high 
latitudes, respectively. The global emission of CH 2 Br 2 is es- 
timated to be 57 Gg Bryr 1 . Note that we derive our emis- 
sion estimates using the background concentration as an ob- 
servational constraint. In general, global chemistry models, 


Atmos. Chem. Phvs., 10, 2269-2286, 2010 


www.atmos-chem-phys .net/ 1 0/2269/20 1 0/ 




Q. Liang et a!.: Finding the missing stratospheric Br y 

even with top-down emission estimate, tend to miss the influ- 
ence of strongly localized sources (with concentrations much 
higher than background concentrations), e.g. along coasts 
and edges of ice sheets. As a result, our global emission esti- 
mate of the oceanic emission for VSL bromocarbons and the 
percentage of coastal emissions (40%) are likely a low limit 
estimate (Quack and Wallace, 2003). The localized sources 
are confined to regions much smaller than a typical model 
grid box (2° latitude by 2.5° longitude in this study). If the 
model were to fill the entire grid box with emissions large 
enough to reproduce the observed localized high concentra- 
tions, the simulated concentrations in the background atmo- 
sphere would be too high compared with observations. 

3.2 Scenario B 

This scenario uses the emission distribution of Scenario 3 
from Warwick et al. (2006). In this scenario, a total of 
380 Gg Bryr~"*for CHBr, and 104GgBryr'' for CH 2 Br 2 
are emitted at the surface. The emissions are concentrated 
in the tropical open oceans, with 75% distributed uniformly 
in the open ocean between 20° S and 20° N and the remain- 
ing 25% between 20° and 50° north and south (Fig. 2 and 
Table 1). 

3.3 Scenario C 

This is based on Scenario 5 of Warwick et al. (2006). We 
follow Warwick et al. (2006) and use a global emission of 
565GgBryr~‘ for CHBr 3 . We reduce the global emission 
of CH 2 Br 2 to 57GgBryr~’ to correct the model high bias 
as compared to observations (Sect. 4). The emission dis- 
tribution contains a combination of open ocean emissions 
and tropical coastline emissions, with ~50% emitted in the 
open ocean between 50° S and 50° N, and the remaining 50% 
along the tropical coastline between 10° S and 10° N (Fig. 2 
and Table 1). 

4 Comparison with observations 
4.1 CHBr 3 

We compare the simulated bromocarbons with the composite 
of aircraft measurements. Figure 4 shows the observed and 
simulated CHBr 3 concentrations in the lower troposphere 
between 800-1000 hPa. We separate aircraft observations 
into two periods, February-May and June-October, which 
approximate the boreal spring and summer seasons, respec- 
tively. Since more than 95% of the observations are made 
between March-May (MAM) and July-September (JAS), we 
compare the spring and summer observations with the sim- 
ulated mean MAM and JAS CHBr 3 concentrations, respec- 
tively. Both the observed and simulated CHBr 3 are regridded 
to a horizontal resolution of 10° x8° for easy comparison. In 
addition, we have added two panels (Fig, 4, middle column) 


2275 

showing simulated CHBr? at locations where there are ob- 
servations available for clear side-by-side comparison. 

High CHBr 3 concentrations are found in a narrow tropi- 
cal band between 10° S-10° N as well as along the coasts, in 
contrast to relatively low values in the northern/ southern cen- 
tral pacific and the center of the continents. This spatial dis- 
tribution is consistent with the short atmospheric lifetime of 
CHBr 3 . As a result, high concentrations of CHBr 3 are found 
only in close proximity to the source regions. CHBr 3 also 
shows significant seasonal contrast in concentrations with 
high values during cold seasons and low values during warm 
seasons. Despite a simple emission scheme with no season- 
ality, the model reproduces well the spatial distribution and 
seasonal variation of CHBr 3 . The simulated CHBr 3 corre- 
lates well with the observations (r = 0.63). This indicates the 
seasonal variability of CHBr 3 is largely controlled by the sea- 
sonality in chemical loss and tropospheric transport (Nielsen 
and Douglass, 2001). The lower concentrations during warm 
seasons are due to efficient chemical destruction by photol- 
ysis and reaction with OH. High emission near the equator 
(between 10° S-10° N) is likely associated with tropical up- 
welling and active planktonic production, as previously ob- 
served in the equatorial Pacific (Atlas et al., 1993) and the 
tropical East Atlantic Ocean (Class et al., 1986). Sensitivity 
simulations show that confining the high open ocean emis- 
sion between 10° S-10° N is essential in re-creating the high 
concentration band in the tropics. This is due to the near- 
zonal advection in the equatorial Pacific associated with sur- 
face trade winds. If the high emission band were extended 
further, an abrupt increase in photolysis rates or OH concen- 
tration from 10° S-10° N to 10-20° S/° N would be needed to 
reproduce the contrast in the observed CHBr 3 mixing ratios 
between 10° S-10° N and 10-20° S/° N, which is unphysical. 
Narrowing the high emission band to 5° S-5° N, on the other 
band, would have led to a high concentration band too nar- 
row compared to the observations when results are plotted at 
2.5° x 2° resolution (not shown). The very high VSLS bro- 
mocarbons observed in the eastern tropical Atlantic (Class et 
al., 1986; Quack et al., 2004 and 2007; Carpenter et al. 2009) 
appear to be a local phenomena. Butler et al. (2007) find 
no evidence for a tropical maxima in CHBr 3 over the open 
Atlantic. We conduct a sensitivity simulation with no tropi- 
cal (between 10° N-10° S) emission maxima in the Atlantic 
to understand how this affects the global emission estimate 
and Br y SLS . The resulted global emission decreases by 1 .6% 
to 418GgBryr~' forCHBr 3 and 55 Gg Bryr^ 1 forCH 2 Br 2 . 
The impact on total bromine in the lower stratosphere is also 
small, ~0.12pptv. 

Figure 5 shows the latitudinal dependence of CHBr 3 at the 
surface. We average aircraft measurements between 0-1 km 
for every 10° latitude band. Surface CHBr 3 maximizes in the 
tropics as well as the high latitudes. Minimum CHBr 3 con- 
centrations are found at ~30° S/° N. The peak in the tropics 
is due to high emissions while the peaks at the high latitudes 
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CHBr 3 


Fig. 4. Left column: the composite of airborne observations of CHBr 3 in the lower troposphere ( 1 000-800 hPa) during boreal spring and 
summer. Middle column: the simulated 1000-800 hPa mean CHBr 3 between March-May (upper panel) and July-September ( lower panel), 
sampled at the same location as the observations. Right column: same as the middle column but for the entire Pacific/North American region. 
The model results are from the simulation run with emission Scenario A. 



Latitude 


Fig. 5. Latitudinal dependence of CHBr 3 at the surface. Black 
line (open squares indicate mean and vertical bar indicate one stan- 
dard deviation) shows aircraft observations averaged below 1 km 
for each 10° latitude band. Surface measurements from Table 2 
in Quack and Wallace (2003) are plotted in gray line with filled 
squares. The simulated zonal averaged CHBr 3 between 0-1 km are 
plotted in color. We use consistent colored lines for model results 
from different simulations for this figure, Figs. 6, 7, 9, and 10: sim- 
ulations with emission Scenario A in red, Scenario B in green, and 
Scenario C in blue. 


are the result of relatively long lifetime. This latitudinal vari- 
ation agrees well with that from previous surface measure- 
ments listed in Table 2 in Quack and Wallace (2003), which 
includes all published atmospheric mixing ratios of CHBr 3 
prior to 2003. The model reproduces relatively well the lat- 
itudinal dependence of CHBr 3 except an underestimate of 
~lpptv in the Arctic. This model-observation difference 
is likely a result of under-representation of emission at the 
ice- water interface, associated with emission from ice algae 
(Sturges et al., 1992, 1997), which is difficult to reproduce 
with a coarse model resolution. 

We further compare the simulated vertical profiles of 
CHBr 3 with observations for individual aircraft missions 
(Fig. 6 ). We sample the model at the nearest grid point 
to the measurement location in the corresponding month. 
Mixing ratio of CHB 13 at the surface is ~ 1-2 pptv, and de- 
creases rapidly to ~ 0.6 pptv in the middle troposphere and 
to <0.4 pptv near 200 hPa. Although regions covered by in- 
dividual aircraft missions differ greatly from each other and 
we deploy a simplified zonal-uniform emission scheme, the 
model in general reproduces reasonably well the observed 
concentration, variability, and vertical gradient during most 
of the missions. The overestimate during ARCTAS-A is 
most likely due to the lack of seasonality in our oceanic emis- 
sion of VSL bromocarbons. To match both the observed pro- 
file from ARCTAS-A (April) and ARCTAS-B (July) would 
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INTEX-B Texas: Mar 



PEMTropics-B: Mar-Apr 



CHBr 3 (pptv) 


TC4: July-August PEMTropics-A: Sep-Oct 



CHBr 3 (pptv) CHBr 3 (pptv) 


Observation 

— Scenario A 

Scenario B 

Scenario C 


Fig. 6. Comparison of the observed and simulated vertical profiles of CHBr3 in the troposphere during individual aircraft missions. The 
upper row shows missions that took place in the northern mid and high latitudes and the lower row shows those in the tropics and southern 
hemisphere. Black lines show the observed mean concentrations with one standard deviation marked by gray shading. The mean CHjBr? 
vertical profiles from three simulations are plotted in color with horizontal bars indicate one standard deviation. Model is sampled at the 
same location as the observations in the corresponding month. 


require a seasonal-varying emission in the high latitudes with 
about zero emission in winter and a peak in summer. This 
seasonality agrees well with several previous studies, sug- 
gesting that CHBr3 released from macroalgae maximizes in 
summer due to enhanced tissue decay and light stimulation 
(e.g. Goodwin et al., 1997; Klick et al. 1993; Carpenter et 
al., 2000). 

The model does not capture the vertical gradient or vari- 
ability of CHBr3 observed during the INTEX-B Texas de- 
ployment. Significant forest fire plumes were sampled dur- 
ing INTEX-B (Singh et al.. 2009). A close examination of 
trace gas measurements along flight tracks shows that high 
CHBr3 mixing ratios are coincident with high levels of trac- 
ers for biomass burning, e.g. carbon monoxide, hydrogen 
cyanide and acetonitrile. This suggests that biomass burning 
is a possible additional source of CHBrr. High mixing ratios 
of CHBr3 associated with African savannah biomass burning 
plumes were documented in Carpenter et al. (2007). The lack 
of biomass burning emissions can also explain the lack of 
variability in simulated CHBr3 mixing ratios compared with 
the high variability observed in the INTEX-A mission, dur- 
ing which significant biomass burning plumes were sampled 
(Turquety et al., 2007; Liang et al., 2007). 


Observed CHBr3 mixing ratios in the UT/LS during recent 
missions, Pre-AVE, AVE, and TC4 (0.2-0. 3 pptv), are sig- 
nificantly higher than that measured during the early STRAT 
mission (<0.1 pptv) (Fig. 7). This is intriguing since long- 
term observations at the surface display little interannual 
variability in CHBr3 concentrations (WMO, 2007). The 
GEOS CCM reproduces well the observed CHB1-3 concen- 
trations in the UT/LS when compared against measurements 
obtained during Pre-AVE, AVE, and TC4 field missions 
(Fig. 7), suggesting a good representation of the atmospheric 
losses and troposphere-to-stratosphere transport of CHB13 
in the UT/LS in the GEOS CCM. The model fails to pro- 
duce the very low concentrations observed during STRAT 
between 80-200 hPa. A similar problem was identified by 
Dvortsov et al. (1999) and Nielsen and Douglass (2001). 
Both attribute this model-observation difference to exces- 
sive convective transport from the boundary layer to the up- 
per troposphere. The fact that GEOS CCM captures well 
the observed profiles during Pre-AVE, AVE, and TC4 im- 
plies that the overestimate during STRAT may not be due 
to excessive convective transport. The difference between 
the CHBr3 during STRAT and that during Pre-AVE has been 
previous noted in Sinnhuber and Folkins (2006) and was sug- 
gested to reflect the differences in the activity of deep con- 
vection between the missions. We have examined in detail 
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Fig. 7. Similar to figure 6 but for the UT/LS region. 


the location of measurements and how the canister samples 
were processed during STRAT. We conclude that the dif- 
ference is likely due to two reasons. First the difference is 
possibly due to the natural variability of short-lived species 
in the upper troposphere and their relationship to recent or 
past convection. The spatial domain of the STRAT measure- 
ments between 80-200 hPa (where there are significant dif- 
ferences between observations and simulated concentrations) 
is much smaller compared to the other missions. Measure- 
ments from the Pre-AVE, AVE, and TC4 missions indicate 
that there is significant spatial variability in the observed con- 
centrations of CHBr 3 , varying from ~0 pptv to >0.4 pptv. 
The STRAT measurements were likely obtained in a rela- 
tively clean pocket of air that was not recently influenced by 
convection and therefore show very low CFlBr 3 concentra- 
tions. Secondly, it is possible that uptake of CHBr 3 on can- 
ister surfaces contributed to lower observed bromoform con- 
centrations in the upper troposphere. Though the canisters 
for STRAT were pre-cleaned with moist air, the treatment 
of sample canisters for the STRAT mission did not include 
addition of extra water vapor that became the standard pro- 
tocol for subsequent missions. The extra water vapor helped 
minimize adsorptive losses to the canister surfaces. 

Due to its longer lifetime, CFEBrj displays less variability 
in its atmospheric concentration compared to CFIBr 3 . This is 
clearly seen in lower tropospheric distribution (Fig. 8) as well 
as vertical profile (Figs. 9 and 10) observed during previous 
aircraft missions. Annual mean mixing ratios of ClDBti are 
~ 1 pptv throughout the troposphere with weak vertical gra- 
dients near the surface (Fig. 9). The mean mixing ratio of 
CFbBri is ~0.5pptv near the tropical tropopause (Fig. 10). 
Most of the variability in the troposphere is tied to the sea- 
sonality in lifetime due to higher rates of reaction with OH 
in warmer seasons. 

The model reproduces well the observed horizontal dis- 
tribution of CH 2 Br 2 in the lower troposphere (Fig. 8 , r = 
0.78) as well as the vertical gradient in the free troposphere 
(Fig. 9). Sensitivity simulations show that, unlike CHBr3, 


CH 2 Br 2 is not sensitive to details in the open ocean vs. 
coastal proportion in emission distribution. This is consistent 
with the fact that CFLBn has a relatively long lifetime and 
is more uniformly mixed. Yet its lifetime is short enough to 
differentiate emissions between tropics, middle latitudes, and 
high latitudes. An accurate latitudinal distribution in emis- 
sion is important in reproducing the atmospheric concentra- 
tion and vertical gradient of GFDB^ at the corresponding 
latitudes. While we did not derive an emission distribution 
specifically for QHB^, our simulated CHiBr 2 and its verti- 
cal gradient matches well with observations. This is consis- 
tent with the fact that both are emitted from the same ocean 
macroalgae and their atmospheric concentrations are highly- 
correlated (Yokouchi et al., 2005). 

Simulated CH 2 Br 2 agrees well with the UT/LS observa- 
tions below 150hPa, but is generally too high compared 
with those between 50-150 hPa, except during Pre-AVE 
(Fig. 10). While the observed CH 2 Br 2 show significant in- 
terannual variation (varying between 0 .3-0.6 pptv at 100 hPa 
from individual missions), the simulated CH 2 Br 2 vary lit- 
tle (<0.1 pptv) from year to year. The model overestimate 
is likely due to the simplification of using the same zonal- 
averaged 2-dimensional monthly mean OH fields for each 
year. However, the fact that simulated CH 2 Br 2 shows a con- 
sistent high bias between 50-150 hPa during most of the mis- 
sions suggests our OH is possibly low compared to mean 
conditions. 

42 Impact of emission 

In this section, we examine the impact of differences in emis- 
sion on atmospheric concentrations of CHBrs and CT^Bn. 

Figure 1 1 compares the annual mean CHBr 3 from year 3 
of simulations A, B, and C. Simulations B and C reproduce 
relatively well the observed high concentration band in the 
tropics, but the band (CHBr 3 > 1 pptv) is too wide compared 
to observations due to reasons discussed above (Sect. 3). The 
simulated concentrations at the mid and high latitudes are 
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Fig. 8. Same as Fig. 4 but for CFbBrs. 
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Fig. 9. Comparison between the observed and simulated vertical profiles of ClLBn in the troposphere. Observations from all available 
aircraft missions are averaged for 30-60° S. 30° S-30° N, 30-60° N, and 60-90° N latitude bands at 1 -km vertical interval. Black lines show 
the mean concentrations with gray shadings indicate one standard deviation. The mean CHbBr 2 vertical profiles from three simulations are 
plotted in color with horizontal bars indicating one standard deviation. Model is sampled at the same location as the observations in the 
corresponding month. 


significantly lower than that from simulation A. This inade- 
quacy in the latitudinal gradient of CHBr 3 is also apparent 
when compared to surface measurements (Fig. 5). Simula- 
tions B and C also have difficulty in reproducing the vertical 
gradient of CHBr 3 for aircraft missions that took place at the 
middle and high latitudes (Fig. 6). These simulation biases 
are due to an underestimate in oceanic emission at the corre- 
sponding latitudes. Sensitivity simulations show that as a re- 
sult of its short lifetime, the influence of surface CHBr 3 emis- 
sions does not extend far, either horizontally or vertically. To 


reproduce the observed vertical profile, it is necessary to have 
oceanic emission underneath or in adjacent regions. The rate 
that CHBr 3 concentration falls off with altitude is sensitive 
to the emission magnitude, and thus is a useful constraint 
to derive emission strength in the corresponding region. In- 
terestingly, while all three simulations have quite different 
emissions, they reproduce well the observed vertical profiles 
of CHBr 3 during PEM-Tropics A and B in the central and 
southern Pacific (Fig. 6). This is consistent with the results 
in Warwick et al. (2006), who found both Scenarios B and C 
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Fig. 10. Same as Fig. 6 but for CH 2 Br 2 in the UT/LS region. 


yield a good match with observations. While the simulated 
CHBr 3 from simulation C agrees well with the tropical open 
ocean observations obtained during PEM-Tropics, it is too 
high compared with the troposphere/lower stratosphere mea- 
surements obtained during TC4, Pre-AVE, and HAVE2 along 
the tropical west coast of South America (Figs. 6 and 7). 
This implies that the 280GgBryr” 1 tropical coastal emis- 
sions used in Scenario C is likely too high. This agrees with 
the results from Kerkweg et al. (2008) who found that with 
emission Scenario 5 from Warwick et al. (2006) (Scenario C 
in the study), the simulated CHBrs is too high compared to 
observations in the UT/LS. 

The simulated CFl 2 Br 2 concentration is sensitive to the 
magnitude of global emission as well as latitudinal distribu- 
tion. Simulation B produces much higher CH 2 Br 2 concen- 
trations compared with the observations and the other two 
simulations in most of the troposphere and LS except in the 
northern high latitudes (Figs. 9 and 10). The global emission 
of 104 Gg Bryr~ 1 used in simulation B was deduced in War- 
wick et al. (2006) by matching the background atmospheric 
concentrations. Kerkweg et al. (2008) also found that their 
simulated CH 2 Br 2 with the global emission estimate sug- 
gested by Warwick et al. (2006) is too high compared against 
the vertical profiles obtained during PEM-Tropics A and B. 
The difference is most likely due to differences in the OH 
fields in individual models. OH in the atmosphere is con- 
trolled by various chemical and moist physical processes and 
an accurate representation of OH in global chemistry remains 
an active yet challenging research task. The difference be- 
tween simulation A and C (Fig. 9) illustrates the sensitivity 
of atmospheric CH 2 Br 2 mixing ratios to emission at the cor- 
responding latitudes. While CH 2 Br 2 displays a much weaker 
vertical gradient compared to CHBrj , its vertical gradient is a 
useful indicator of surface emissions, similar to CHBr 3 . With 
little emissions in the mid and high latitudes (simulation C), 
it is difficult to reproduce the observed concentrations as well 
as the vertical gradients. 


The above suggests that top-down emission estimate of 
short-lived bromocarbons are subject to great uncertainties 
introduced by the abundance of available measurements, 
both spatially and temporally, and model OH fields. For 
CHBi' 3 , which has a relatively short lifetime and loss occurs 
dominantly via photolysis, the derived estimate is depen- 
dent on the observations used. The uncertainty can be sig- 
nificantly reduced when there are more measurements avail- 
able, particularly over the strong emission regions, e.g. along 
the coasts and near the equatorial tropics. Dibromomethane 
has a relatively long lifetime and is well-mixed in the tro- 
posphere, with loss occur predominantly via reaction with 
OH. Thus, uncertainty in the top-down emission estimate 
of CH 2 Br 2 is mostly associated with the calculated rate of 
loss against OH. Differences between this study and that 
from Warwick et al. (2006) suggest that the uncertainty can 
be as high as 100%, most likely due to differences in OH. 
Top-down emission estimates for CHBrCL, CHBr 2 Cl and 
CH 2 BrCl are subject to similar uncertainties as their loss oc- 
curs dominantly through reaction with OH. Results from this 
study also indicates that while the observed vertical gradient 
is a useful constraint in deriving surface emissions as shown 
by many previous studies, it is helpful to combine vertical 
profiles with horizontal distributions to achieve an accurate 
emission estimate globally and a better understanding of pos- 
sible model bias. 


5 The contribution of short-lived bromocarbons to 
stratospheric bromine 

The simulated CHBr 3 and CH 2 Br 2 concentrations in the 
troposphere reach steady state in about or less than one 
year. In contrast, their impact on total bromine levels in the 
stratosphere, particularly Br y , approaches steady state much 
more slowly, a result of slow transport in the stratosphere. 
Figure 12 shows the temporal increment of total bromine 
(organic and inorganic) produced from VSL bromocarbons 
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(a) Emission Scenario A 



(b) Emission Scenario B 
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Fig. 11. The annual mean horizontal distribution of CHBr 3 in 
the lower troposphere ( 1 000-800 hPa) from simulations conducted 
with emission (a) Scenario A, (b) Scenario B, and (c) Scenario C. 

(Br VSI -S) at different locations in the stratosphere. Br VSLS 
near the tropical tropopause (0°N, lOOhPa) reaches steady 
state concentration at the end of year 2 while that near the 
mid/high latitude tropopause and in the tropical mid/upper 
stratosphere gradually increases. The entire stratosphere ap- 
proaches steady .state at year 12. The zonal averaged Br VSLS 
from year 3 of the simulation shows that the troposphere-to- 
stratosphere transport of Br VSLS occurs via vertical advection 
into the tropical stratosphere and quasi-horizontal isentropic 
transport into the mid-latitude lower stratosphere (Fig. 13c). 
similar to results from Nielsen and Douglass (2001). The 
middle and high latitude lower stratosphere is controlled by 
the relatively rapid quasi-horizontal cross-tropopause isen- 
tropic transport from the tropical troposphere and the slower 
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Fig. 12. Simulated total bromine at 100 hPa (gray solid line). 50hPa 
(gray dashed line), and 10 hPa (gray dash-dotted line) at the equator, 
as well as those at 45° N (red) and 90° N (blue) at 100 hpa between 
January 1990 and December 2001. We apply a 13-month running 
mean filter to each monthly-mean time series to remove seasonal 
variations. 


descending motion associated with the mean meridional cir- 
culation (Nielsen and Douglass, 2001). As a result, Br VSLS in 
the middle/high latitude lower stratosphere increases rapidly 
at the beginning of the simulation and approaches steady 
state more slowly than that at the equatorial middle and high 
stratosphere (Fig. 12). 

Figures 1 3d— f and 14 present the steady-state contribution 
of CHBrs and CFLBn to bromine in the atmosphere. At 
steady state, inclusion of CHBr 3 and CFLB^ adds a uni- 
form ~5pptv to total bromine in the entire stratosphere, 
set by Br VSLS concentration at the entry point - the tropi- 
cal tropopause layer. Previous modeling studies with 3-4 
years of simulation length suggested that contribution of very 
short-lived bromocarbons to inorganic bromine in the strato- 
sphere peaks at a certain altitude in the stratosphere (e.g. 
Nielsen and Douglass, 2001: Warwick et al., 2006). How- 
ever, our results demonstrate that this is merely a snap-shot 
contribution amid the spin-up process. As we can clearly see 
in Fig. 13b, at year 3, Br y SLS peaks in the tropical lower 
stratosphere. As the simulation continues, the maximum 
BrVSLS area extenc i s further both horizontally and vertically 
until reaching steady state concentration throughout the en- 
tire stratosphere. 

Troposphere-to-stratosphere transport of source gases and 
their degradation products are equally important. Inorganic 
bromine accounts for half of Br VSLS near the tropopause 
and its contribution increases rapidly with altitude as source 
gases quickly photo-dissociate to produce Br y . Above 
50hPa, almost all is in the form of Br y . The contribution 
of CH?Br 2 to reactive bromine in the stratosphere is of equal 
or greater importance as that of CHBr 3 . While CHBr 3 con- 
tributes mostly via product gas injection, source gas injection 
is a more important pathway for ClLBr;;. This additional 
Br y produced from the VSLS can have a significant impact 
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Fig. 13. The latitude-pressure cross section of the contribution of CHBr 3 and CH 2 Br 2 to organic bromine (left column), inorganic bromine 
(middle column) and total bromine (right column) from year 3 and year 1 2 of the simulation. 


on stratospheric ozone depletion. The associated BrO cause 
more ozone depletion due to enhanced catalytic loss by the 
BrO+CIO cycle (Salawitch et al., 2005). Secondly, an ad- 
dition of ~5 pptv of Br y causes the photochemical loss of 
ozone by BrO+H 02 cycle to be of considerable importance 
below ~14km (Salawitch et al., 2005). 

About 30% of the Br y produced from CHBr 3 and CHtB^ 
in the troposphere is removed by wet scavenging, and more 
than 85% of the wet removal is due to large-scale pre- 
cipitation below 500 hPa. The mixing ratio of Br VSLS 
is relatively constant from the tropical middle troposphere 
to the stratopause, implying small wet scavenging during 
troposphere-to-stratosphere transport. We conducted a 12- 
year sensitivity simulation without wet scavenging in con- 
vective updrafts and found that convective scavenging only 
accounts for ~0.2pptv (4%) difference in BrJ SLS in the 
stratosphere, contrary to the conventional wisdom. This 
result agrees with the sensitivity study of Dvortsov et 
al. (1999), who found that estimates of Br y in the lower- 
most stratosphere are not very sensitive to convection. Mix- 
ing ratios of a chemical species entering the stratosphere 
is determined by its photochemistry and dehydration pro- 
cesses during troposphere-to-stratosphere transport through 
the TTL (Sinnhuber and Folkins, 2006; Fueglistaler et al., 
2009). Based on the tropical vertical profiles of Br^ SLS 
(Fig. 14), we can deduce that about 3 pptv of the strato- 
spheric Bry SLS is formed either in the stratosphere or in 
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the TTL above the level of neutral buoyancy (LNB) (poten- 
tial temperature ~355K and ~125hPa) where they can en- 
ter the stratosphere via radiative ascent (Fueglistaler et al., 
2009), and are therefore not subject to convective scaveng- 
ing. The remaining ~2pptv are formed below the TTL. 
The results from our detailed wet deposition scheme, which 
includes rainout/washout, convective scavenging and evap- 
oration, suggests that dehydration during troposphere-to- 
stratosphere transport occurs slowly so that most of the 
BrVSLS j s i 0 f tet | j n convective updrafts is released back 
to the atmosphere during evaporation, subsequently escaping 
to the stratosphere. 

Our estimate of ~5 pptv contribution of CHBr 3 and 
CFLBr 2 to stratospheric bromine agrees well with previous 
ground and balloon-borne measurements (~5 pptv) (Sinnhu- 
ber et al., 2002; Dorf et al., 2008) as well as satellite observa- 
tions (3-8 pptv) (Sinnhuber et al., 2005; Livesey et al., 2006; 
Sioris et al., 2006). This estimate is similar to that from War- 
wick et al. (2006) (6-7 pptv for all five VSL bromocarbons), 
but in general higher than many previous modeling estimates 
that between 0.7-1 .8 pptv for CHBn alone (Dvortsov et al., 
1999; Nielsen and Douglass, 2001 ; Hossaini et al., 2010) and 
2.4-2.8 pptv for CHBr 3 and CFbBr 2 combined (Gettelman et 
al., 2009: Hossaini et al., 2010). The difference is mostly due 
to differences in the rate of wet scavenging, which is ~50% 
lower in the GEOS CCM (wet deposition lifetime ~15 days) 
than most models (wet deposition lifetime of 10 days). While 
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Fig. 14. The model calculated contribution of organic bromine (red solid lines), inorganic bromine (blue lines), and total bromine (black 
lines) from short-lived bromocarbons, averaged between (a) 30° S-30° N, (b) 30-60° S/N. and (c) 60-90° S/N. The relative contribution of 
bromine from CHBr 3 (CHBrj x3, red dotted lines) and CFLB^ (CILBr? x2, red dashed lines) to organic bromine is also included. 


Warwick et al. (2006) uses the same wet deposition rate as 
the other models, the impact of higher emissions cancels off 
with the impact of faster scavenging and yields a comparable 
contribution as this study. Similarly, our estimated contribu- 
tion due to injection of product gases which are subject to wet 
scavenging (~50% for CHBr 3 and CHLBra combined) is also 
significantly higher than the ~50% for CHBr 3 (Dvortsov et 
al., 1999; Nielsen and Douglass, 2001 ; Hossaini et al., 2010) 
and ~90% for CFLB^ (Hossaini et al., 2010) from previous 
modeling analysis. 


6 Conclusions 

We performed a modeling study using the GEOS CCM to 
quantify the contribution of CHBr 3 and CH 2 Br 2 to reactive 
bromine in the stratosphere. We conducted a 12-year bro- 
mocarbon simulation that includes a detailed representation 
of oceanic emissions of CHBr 3 and their chemical 

losses through photolysis and reaction with OH, surface dry 
deposition of inorganic bromine (their degradation product), 
and wet scavenging of inorganic bromine in large-scale pre- 
cipitation and convective updrafts. 

We used observed concentrations and vertical profiles 
of CHBr 3 and Cl-LBt^ from previous NASA aircraft cam- 
paigns, including PEM-tropics, TRACE-P, INTEX, ARC- 
TAS, TC4, STRAT, Pre-AVE and AVE, to derive a top- 
down emission estimate of CHBrj and CHtB^ from oceanic 
sources. To produce the observed background mixing ra- 
tios, we estimated that the global emissions of CHBr 3 and 
CTLBd are 425GgBryr -1 and 57 GgBryr -1 , respectively, 
with 60% from open ocean and 40% from coastlines. Our 
emission estimates agree well with previous modeling es- 
timates and are within the uncertainty range of bottom-up 
emission estimates. The 425 GgBryr -1 for CHBr 3 is likely 
a lower limit as the model does not reproduce the very high 
mixing ratios observed in the marine boundary layer. Our 


sensitivity simulations indicate that top-down emission esti- 
mate for CHiBn is sensitive to model representation of at- 
mospheric OH concentration as reaction with OH is the dom- 
inant process for degradation of CH 2 Br 2 . 

When initialized with zero atmospheric concentrations, 
the distribution of bromocarbon source gases and the re- 
sulted inorganic bromine to the stratosphere approach steady 
state in ~12 years. At steady state, including CHBr 3 and 
CH 2 Br 2 adds a ubiquitous ~5pptv to total bromine (or- 
ganic and inorganic) in the stratosphere, with transport of 
source gases and transport of their degradation product be- 
ing equally important. Inorganic bromine accounts for a half 
(~2.5pptv) of total bromine from the VSL bromocarbons 
near the tropopause. Its contribution increases rapidly with 
altitude and reaches a constant ~5 pptv at 50 hPa and above, 
which can have a significant impact on stratospheric ozone 
depletion. 

Bromoform has been suggested to play the most important 
role among short-lived oceanic bromocarbons in contributing 
reactive bromine to the stratosphere. Our results indicate that 
CHiBr 2 is of equal or greater importance. While CHBr 3 con- 
tributes mostly via product gas injection, source gas injection 
is a more important pathway for CHiBr 2 . 

More than 85% of the wet scavenging of Br^ SLS occurs 
below 500 hPa due to large-scale precipitation. Great em- 
phasis has been placed on the impact of convective trans- 
port on the troposphere-to-stratosphere transport of the VSL 
bromocarbons and wet scavenging of the resulted inorganic 
bromine. However, our sensitivity study shows that con- 
vective scavenging only accounts for ~0.2pptv (4%) differ- 
ence in inorganic bromine delivered to the stratosphere due to 
dominant source gas injection above neutral buoyancy level 
in the tropical tropopause layer and slow dehydration during 
lofting in deep convection. 
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